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Abstract. We study a class of slow-fast Hamiltonian systems with any finite number 
of degrees of freedom, but with at least one slow one and two fast ones. At e = the 
slow dynamics is frozen. We assume that the frozen system (i.e. the unperturbed fast 
dynamics) has families of hyperbolic periodic orbits with transversal heteroclinics. 

For each periodic orbit we define an action J. This action may be viewed as an 
action Hamiltonian (in the slow variables). It has been shown in [4] that there are 
orbits of the full dynamics which shadow any finite combination of forward orbits of 
J for a time t = 0{£~'^). 

We introduce an assumption on the mutual relationship between the actions J. 
This assumption enables us to shadow any continuous curve (of arbitrary length) in 
the slow phase space for any time. The slow dynamics shadows the curve as a purely 
geometrical object, thus the time on the slow dynamics has to be reparameterised. 
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1. Introduction 

Let us consider the Hamiltonian system defined by the Hamiltonian function 

H = H{x,y,u,v;e) 
and the symplectic form 

Q = dy A dx -\ — dv A du 

e 

where x,y & M^™" and u,v & M^''. We assume that the parameter e is small, hence the 
Hamiltonian system is "slow-fast" with {x,y) being fast variables and {u,v) being slow 
ones. This can readily be seen from the equations of motion 

(1) 

We note that the form of (1) is not unusual, it appears in many applications. 
If we set £ = the equations of motion are reduced to 



• _ aff M = 0, 



(2) 



We refer to (2) as the frozen system. 

We are interested in describing the slow dynamics of system (1). In general 
it is a difficult task to exactly describe the slow dynamics and therefore often only 
approximations, or averaged solutions, arc sought. The procedure to obtain the averaged 
solutions, i.e. an averaging method, is typically tailored for the particular class of 
systems it is applied to. The averaging method for the classical problem where there is 
only one fast degree of freedom is described in [3]. This method has been generalised 
to systems with several fast degrees of freedom where the fast system rotates with a 
constant vector of frequencies, see [10]. There are also averaging methods for slow- fast 
systems where the fast system is uniformly hyperbolic, see [1] (or [8] for a description in 
English). Contrary to [10] the class of problems studied in [1] contains systems which 
are fully coupled, i.e. the fast system depends on the slow variables as well as the fast 
ones. Or, more generally, if the fast system is ergodic and satisfies assumptions on how 
fast time averages converge to space averages then, see [7], the fast dynamics can be 
averaged out. Regardless of which averaging method is employed the aim is to derive 
effective equations for the slow dynamics which arc independent of the fast variables 
{x, y) . We may write such averaged slow dynamics as 

/dH\ . /dH\ 

where (■) denotes averaging with respect to {x,y). A crucial point in the justification 
of all averaging methods is to verify that time averages can be approximated by space 



Geometric shadowing in slow-fast Hamiltonian systems 



3 



averages. This property holds for most but not all trajectories. For example, if the fast 
dynamics spends most of its time close to a periodic orbit there is, in general, no reason 
to believe that the time average would converge to the space average. 

In fact when the fast dynamics spend most of its time close to periodic orbits the 
slow dynamics of (1) may behave very differently from the averaged dynamics (3). In [4] 
it was shown that if the system satisfies assumptions [Al] and [A2] (defined later) then 
there exist trajectories of the full dynamics whose slow component deviates significantly 
from (3). In particular it was shown that there are orbits shadowing accessible paths of 
finite length composed of forward trajectories of the auxiliary systems 

dJc . dJc 

ov ou 

where Jc is an action defined on a periodic orbit (labelled by c) in the fast phase space 
(sec Definition 1). This shadowing result is valid on time-scales of order 0{£~^). We 
note that this method yields trajectories which deviate at the rate 0{£) from (3). 

This approach is a generalisation of the mechanism proposed in [5] for studying drift 
of the energy in a Hamiltonian system which depends on time explicitly and slowly. In 
this set-up it was shown that switching between fast periodic orbits does indeed provide 
the fastest possible rate of energy growth in several situations (see [5]). An interesting 
direct application of this theory is a rigorous proof of Fermi acceleration for a class of 
billiards with slowly moving boundary, [6]. 

The main result of this paper eliminates the upper bound on the time for which the 
shadowing result holds, that is we provide a description of the shadowing orbits for all 
times. Moreover, we show that for any continuous curve in the slow phase space there is 
a trajectory of the full dynamics whose slow component shadows it. We achieve this by 
refining the mechanism in [5] and [4] such that for any two 0{e) close points in the slow 
phase space there is a trajectory which starts in a neighbourhood of the first point and 
ends up in a neighbourhood of the second one. To prove this we have to introduce an 
assumption on the mutual relationship between the actions Jc- This assumption ensures 
that any vector in the tangent space of a point in the slow phase space can be written 
as a linear combination of gradients of the actions Jc where all coefficients are positive. 
We identify this linear combination with a path in the slow phase space and denote it 
guiding path (which is a generalisation of accessible path). We only consider guiding 
paths of length 0{e) and show that for any such guiding path there exists a trajectory 
of the full dynamics whose slow component shadows it. Thus by approximating any 
given curve by a set of points which are 0{e) apart we can find trajectories of the full 
dynamics whose slow component shadows a guiding path between these points. Then 
[5] implies that that there exists a trajectory which lies close to the union of trajectories 
shadowing the guiding paths, hence there is a trajectory which shadows the entire curve. 
The idea of the proof is similar to that of [4] , that is, using that the full dynamics takes 
place on normally hyperbolic manifolds where certain action Hamiltonians are preserved 
for long times. 
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The paper is structured as follows. In the next section we present the assumptions 
on the frozen system (2) and state our main result. Before proving this result we present 
an example where the assumptions on the frozen system are verified to hold. In Section 
3 we provide a short summary of the results from [4] and [5] on normal hyperbolicity 
which we will require to prove our main theorem. In Section 4 we state two lemmas 
which we then combine to prove Theorem 1. The proofs of the two lemmas are of more 
technical nature and have been postponed to Section 5 to increase the readability of the 
previous section. 

2. Set-up and Statement of the Result 

We impose a number of assumptions on the frozen system (2). Let D C M^'^ be an open 
and bounded subset. 

[Al] We assume that the frozen system (2) has n families of hyperbolic periodic orbits 
Lc{u, v), c e {ci, . . . , c„} , defined for all v) G D. 

[A2] We assume that each of the periodic orbits has a family of hetcroclinic orbits to 
every other periodic orbit, i.e. for all Cj, Cj € {ci, . . . , c„} and [u, v) & D there are 
a pair of transversal heteroclinic orbits 

T,^,^{u,v) C W^{L,Xu,v))nW'{L,^{u,v)), 

r,,cAu,v) c w^{L,.{u,v))nw\L,,{u,v)). 

We note that under assumptions [Al] and [A2] the frozen system has a family of 
uniformly hyperbolic invariant transitive sets A(u ,„) , also known as Smale horseshoes (see 
for instance [9]). The dynamics on the Smale horseshoe can be described by symbolic 
dynamics. Let A := ^(u,v)eD-^{u,v)- For each family of periodic orbits we define a family 
of actions. 

Definition 1 The action Jc of a periodic orbit Lc is defined by the integral 
Jc{u, v) := (j) ydx. 

J Lc{u,v) 

The function Jc(u, v) is independent of the fast variables and can be considered as 
a Hamiltonian function which generates some dynamics in the slow variables 

• - 1 9Jc . _ 1 dJc , . 

~ T^(u, v) dv' " T,{u, v) du ' 

where Tc is the period of the periodic orbit Lc- Following the notation of [11] we refer 
to (4) as the guiding system and its equations of motion can be written in the concise 
form 



where Xc is the guiding vector field 



(5) 
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[A3] We assume that there is a closed subset V C D C M?'^ such that for all {u, v) e T>, 
is inside the convex envelope (or convex hull) of {V(„,^) Jc-} where i — 1, . . . ,n. 

Reiiicirk 1 Although n — 2d + 1 actions are sufficient to satisfy [A3] we allow for a 
larger set of actions to be considered since it may increase the size of the domain V. 
Nevertheless, at any point in T> we will use only 2d-\- 1 vectors (but not necessarily the 
same 2d-\- 1 vectors for every point in T>). 

We can now state the main result. Let tt : ]R2m+2(i _^ -^d j-^g ^j^g projection on the 
slow variables. 

Theorem 1 Let 7 : M"*" V he any continuous curve. Assume that the frozen system 
(2) satisfies assumptions [Al],[A2j and [A3]. Then there exist positive constants C and 
£0 such that for any < e < £0 there exists a trajectory ^l){t) of the full system (1) and 
a continuous monotone reparameterisation of time T{t) such that the slow component 
z{t) :— 7rip{t) of ip{t) satisfies z{0) — 7(0) and 

\\z{T{t))-j{t)\\<Ce, 

for all t > 0. 

Before proving the theorem we construct an example where assumptions [Al], [A2] 
and [A3] are shown to be satisfied. 

2.1. Example 

Let £ = and let be a small parameter. Consider a Hamiltonian function of the form 

H{x, y) = Ho{x, y) - iiHi{x, y) 

where (x, y) e M^. Assume that the dynamics of the "unperturbed" Hamiltonian 
has three hyperbolic periodic orbits Lc, c e {ci, 02,03} with transversal heteroclinic 
connections (hence assumption [Al] and [A2] are satisfied). The implicit function 
theorem implies that the hyperbolic periodic orbits persists for jj, sufficiently small, and 
furthermore the periodic orbits depend smoothly on Therefore we can expand the 
periodic orbit as well as the action Jc in a power series in //. For the action we write 

J,^r^ + ^iJl + o{^^), 

where Jq is the action of the periodic orbits of the "unperturbed" Hamiltonian Hq. Jq 
and Ji are constants as the Hamiltonian is independent of {u,v). The action Jq and its 
first order correction Ji are generically non-zero. Now, by abusing the notation slightly 
we let fj, depend on {u, v) in the following way 

3 

H = n{x,y,u,v) = i2^XLc^+s{x,y)^^'{u,v) 

i=l 

3 

+ I^^x{l,^+s) {x,y)f{x,y), 

i=l 
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where is a small parameter, x is the indicator function, Lc- + 5 is a 5 neighbourhood 
of Lc., i.e. 

Lci + 5 := {{x,y) s.t. dist{{x,y), L^i) < 6} 

and / is chosen as to interpolate Hi into a C°° function. Then the action of the periodic 
orbit Lc becomes 

Jciu,v) = + iiip^'{u,v)Jl' + 0(^2). 

For example, by choosing 

Lp"^ = sgn{Jl^)v, 

ip"^ = - sgn{Jl^)u, 

(^■=3 — sgn{Jl^)u — sgn{Ji^)v, 

where sgn{A) := 1 if ^4 > and sgn{A) := — 1 if ^4 < 0, will generate 

This choice ensures that assumption [A3] is satisfied. 
3. Normal Hyperbolicity and symbolic dynamics 

It was shown in [4] and [5] that the dynamics of the full system (1) can be studied using 
symbolic dynamics. Here we provide a short summary of the results therein which we 
will need to prove Theorem 1. For each periodic orbit Lc we denote by Ec a Poincare 
section, and by Xi,yi,Zi (where Zi = {ui,Vi)) we denote an intersection of a solution of 
(1) with one of these Poincare sections. Let ^ = {^iy~°°^ be a bi-infinite sequence of 
letters G {ci, . . . , c„} where q is the index of the n actions. The sequence ^ is called 
the code of the dynamics. It has been shown that by specifying an initial condition for 
the dynamics the code ^ generates a trajectory of (1) whose dynamics on the Poincare 
sections are given by 

Zi+i = Zi + e(f)^^^^^^{xi{zi,^,e),yi+i{zi,^,e), Zi,e), (6) 

where G C^, and the index i counts intersections with the Poincare sections. The 
iterate Zi G S^.. Note that the functions Xi and yi depends on entire code ^. 

Lemma 2 in [4] implies that for £ = and any two codes ^^^^ and ^^"^^ that satisfy 
^j^^-* = for \i\ < n for any n the following estimate holds 

me.x{\\x.{z,e'^)-x.{z,^^'^)\\ , Wmiz,^'^) - y,{z,C^'^)\\} < 2rA"-l^ (7) 

where the constants r > and < A < 1 do not depend on the sequences 

^(i),(2). For 

e small Lemma 3 in [4] implies the functions {xi,yi) are defined for all small e and all 
z & V, they are uniformly bounded along with their first derivatives with respect to z 
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and satisfy (7). Moreover, by the lemma, there is a constant Co > 0, independent of 
the code ^, such that 

\\xi{z,i,e)-Xi{z,i,Q),yi{z,i,e)-yi{z,i,Q)\\ < CqE, (8) 

for all i e Z. 



4. Proof of Theorem 1 

In this section we construct a code ^^*^ that generates a trajectory z^*^ of the full 
dynamics (6) which satisfies Theorem 1. The proof relies on two lemmas. The first one 
states that two trajectories of (6) which have a slow iterate in common stay uniformly 
close to each other as long as their respective codes coincide (Lemma 1). The second one 
states that given any two points in the slow phase space which are 0(e) close we can find 
a trajectory which goes from the neighbourhood of the first point to the neighbourhood 
of the second point. This trajectory is essentially obtained by updating the code of 
another trajectory. We then combine these two lemmas to give an inductive proof of 
Theorem 1. 



4.1. Uniform closeness of trajectories 

As mentioned in the previous paragraph, finding a trajectory that satisfies Theorem 1 
is an iterative process which involves updating the code. The crucial point here is that 
when we update the code we do not only alter the future of the trajectory but we switch 
to another trajectory. Differently put, as we update the code the entire slow dynamics 
changes (not only the "future" iterates of the trajectory); in fact if one considers the 
description of the slow component of the full dynamics given in equation (6) one sees 
that the first two arguments of the function 4'e,i£.i+i depend on the entire code ^ (not 
just the current and next to current code elements). The following lemma gives us a 
uniform estimate on how much the "past" of the trajectory changes when the code is 
updated and the initial condition is kept fixed, i.e. z^^^ = z^q\ The estimate appears 
in [5] (see the paragraph between equation (56) and (57) therein) but since the proof 
was only sketched and the result is crucial to our theory we state it together with a full 
proof. 

Lemma 1 There exists £0 > such that for any Q < e < sq the following holds: If ^^""^ 
and ^^''^ are two codes such that for some N > 

and z^""^ and z^''^ are two slow trajectories of (6) which correspond to the codes and 
^(h) respectively and z^^ = Zq^ then 

< Ke for all Q<i<N (9) 



Ja) _ (b) 



where 



K:=^^. (10) 
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Note that the constant K is independent of N. We postpone the proof of this 
Lemma to Section 5.1. 

With K given by (10) we let 

4r 

Ai = max sup 3£(X+ ||(/)cc||c;i t + + \Tc{u,v)\\Xc{u,v)\) 



and define recursively 



Ai :— max sup 3e{Ai_i + ||0, 
^ {u,v)ev 



Ar 



cc\\c^ 



A 



+ l + \T,{u,v)\\X,{u,v)\). 



Let 



where is the number of slow degrees of freedom. 



(11) 



Lemma 2 Let L > be any fixed constant and A given by (11). There exists eo > 
such that for any < e < £o CLnd any z E V and any code with a corresponding 
trajectory z^f' satisfying 



< e{L + A) 



the following holds: There exists another code ^^^^ with 

ef^ = ef' ^<o 

such that for any p < there exists N E N and a corresponding trajectory z^''^ which 
satisfies 

P ~ p ' 



< sA. 

Moreover there exists a uniformly bounded constant Ci > such that 

<e{Ci + A) forO<i<N. 



z) ' - z 



(12) 



Next we combine these two lemmas to show that we can shadow any curve 
7 : M+ ^ P. 

4-2. Combining the results 



Proof, (of Theorem 1) Let 7 : — > "D be any curve. Pick a constant L > 
arbitrarily. Take the smallest £0 of Lemma 1 and Lemma 2. The following analysis is 
valid for any < e < £o- Using L we define a sequence tj, z e N, as follows 

(13) 

ti+,^min{t:h{t)-^{ti)\\^eL} 
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and if for some k 

\\l{t) - l{tk)\\ < sL foralH>tfe 

then 

ti^i = ti + 1 for all i > k. 

The sequence ti divide the curve 7 into points 7(^1) which are at most sL apart. Next 
we take any code and a corresponding trajectory z^^^ such that ^q^'' = 7(^0) — 7(0)- 
We also define P(0) = 0. 

Inductive assumption: There exists a code ^^''^ and a monotone sequence P{1) 
such that for all 1 < / < A; 



= er'' forj<p(/-i). 

Moreover, there is a trajectory z^^^ corresponding to ^^^^ such that 



,(/-!) 



.(0 



.(0) 



and 



Furthermore 



,(0 



< sA. 



(14) 
(15) 



114 - 7(^011 < + A) for Pil -l)<j< P{1), 

where Ci is given by Lemma 2. 

Inductive step. We will use Lemma 2 to verify the inductive assumption. Set 
p — —P{k), z — ^(tk+i) and — ^plk)+i ^^^^ using (14) and (13) we get 



< e{A + L). 



Then applying Lemma 2 implies that there exists a code Cp(^)^|j ^i"' which satisfies 

^(^+1) = ^(fc) for i < P{k) 
and generates a trajectory z"^^^^ which satisfies 

^^^^ = ^^=7(0). 
Furthermore there exists A^jt+i such that 



(fc+i) 

P(fc)+Affe+i 



< eA 



and 



fc+1 



l{tk+i)\\ <e{C + A) ior P{k)<j <P{k + l). 



The monotone sequence P{k) is defined inductively by 
P(0) = 0, 

P{k) - P{k - 1) + Nk. 
This concludes the inductive step. 
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By induction there exists a unique code ^(*) which for all > satisfies 
Denote by z^*^ the trajectory of (6) which corresponds to the code ^^*^ and satisfies 



.(*) 



7(0). Next we want to apply Lemma 1 to show that for all k the trajectory z 



(*) 



hes close to the points ^{tk). Set = ^^^^ = i^^^ and N = P{k). Then Lemma 1 
implies 



Z: 



{*) _ Jk-) 



<eK for < 2 < P{k). 



Combining this estimate with (15) gives 



V/O 



4*^ - l{tk) < e{K + Ci + A) for P{k-l)<i< P{k). (16) 



To conclude the proof we have to define the time reparameterisation T{t) of 
Theorem 1 and pass from the discrete solution zf^ to the time continuous z^*\t). Let 
us begin with the time reparameterisation. For <t < i^+i we define 



Tk{t) = Tp(k) + 



Tp(fe+1) — Tp(fc) 



{t-tk), 



tk+1 — tk 

where is the time z^*\t) intersects the Poincare surface E (*>, i.e. z^*\Ti) — zl*\ The 
complete T(t) is obtained by gluing together all Tk(t). 

Next for tk < t < tk+i, and consequently P{k) < i < P{k + 1), we estimate the 
distance between z''*\t) and ^{t) as follows 

||^W(r(t))-7(t)|| 



< ||zW(r(t))-zW(r(4^^))||+ .w(r(4+o)-4*(U 



(*), 



+ 



%+i)-litk+i) +||7(^fe+i)-7WII- 



Consider the right hand side. The second term is by definition, the third term is 
bounded by (16) and the fourth one by (13). The first term in the right hand side we 
estimate as 

\\z^*\T{t)) - z^*\r{tk+i))\\ < sup_isup(rp(fe+i) - Tp(fe)). 

We note that the set of periodic orbits and transversal heteroclinics is compact. The 
full trajectory, whose projection on the slow phase space is z^*\ lies in a compact 
neighbourhood of this set, therefore z given by (1) is uniformly bounded on this set: 
there exists a constant C2 > 0, independent of k, such that 

\\z^*\r{t))-z^*\r{tk+i))\\<eC2. 

By choosing C := C2 + L + X + Ci + ^4 we have 
\\z^*\T{t)) - 7(i)|| < eC for t > 0. 
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In this section we include the proofs of Lemma 1 and Lemma 2 that we used to prove 
Theorem 1 



5.1. Proof of Lemma 1 



Let Eq be sufficiently small for the normal hyperbolicity estimates to be valid and smaller 

l-A 
2CA' 



than ^T=7T, where 



C = sup ( 1 + max 

{u,v)eD 



dx 



dz 



dy 



dz 



(17) 



Consider the two codes and By assumption we have 



_ c(&) 



(18) 



\i\ < N, 

and as a consequence of normal hyperbolicity that 

||x.(z,e('^)) - x,(z,e(')),y,(z,e(»)) - miz,^^'y)\\ < 2rA^^'-l'l 

for all \i\ < N, see (7). For convenience we repeat equation (6) which gives the slow 
component of the full dynamics generated by the code ^ 

Zi+i = Zi + e(f)^.^.^^{xi{zi,e;^),yi+i{zi,e;^),Zi,e). 

We note that the functions Xi and yi depend on the full code ^ (not only the current 
element in the code). Taking the difference of the z components of the two trajectories 
coded by and ^^'^^ and using that they have the same initial condition ^q"^ — Zq^ we 
get 



i=0 



,.u^.>{x,{zf\e;e\yi+,{zf\e-,e\zf\e 



ip) 



Now, i^f' = ^l"' for all \i\ < N, therefore (f) ja)ja) and (j) (b)Jh) are the same functions 
in the interval we are studying. Using the mean value inequality and suppressing the 
dependence of e in the notation we get 

Y,(x,{zt\e<^))-x,{zf^^e\ 



-k ^k 



< e 



KH-.(^l"',e<"*)-K+i(4",«'": 



+ 



i=0 
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By adding and subtracting terms and using the triangle inequality we rewrite this 
expression to a form where inequality (18) can be used 



fc-i 



i=0 



+ 



i=0 

k-l 



+ 



i=0 



Using inequality (18) to estimate the first term and the mean value inequality to estimate 
the second term gives us 

fc-i fc-i 



Ja) _ (b) 
'^k ^fc 



< 



e\\<P\yJ2'^rX''~' + '\\^\yIl 



^1 ^1 



i=0 



i=0 



max 



dx 



dz 



dy 



dz 



1=0 



Aa) _ (b) 



which we rewrite as 

Sa) (b) 



^k ^k 



\N-k+l 

i=0 



where C is given by (17). Now, in order to prove Lemma 1 we will show that 



^k ^k 



(19) 



(20) 



for all < A; < TV where K is given by (9) . The Gronwall type of estimate follows from 
inequality (19) by finite induction. Since z^'' — Zq^ the statement is valid for k — 0. Let 
us assume that (20) is true for all k < m. Consider the case = m + 1. By (19) we get 



<2eU\\c,rr^ + eCY, 

i=Q 



L(a) _ (b) 



The induction assumption implies that 



\N-m 

<2eU\yrj— + sCj2^KX 

i=0 

\N-m \N-m 



N-i 



< 2£||0||^ir^_^ 



1- A 



For 



_ 8Mg.rA 
1 - A 
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and < e < we have 



K > 



rX 



l-X-eC\ 



Thus it follows that 



(21) 



Since < A < 1 and < m + 1 < A'" the lemma follows from (21) . We note that K is 
independent of N and e. □ 



5.2. Proof of Lemma 2 

Fix L > arbitrarily and let A be given by (11). Then we choose Eq as the smallest of 
that required by Lemma 1 and 

1 



where C is given by (17). Then for any e < Sq let ^^""^ be a code generating a trajectory 
of (6) such that 



(a) 

zy -z 



< e{L + A). 

To show that there exists a code ^^^^ that satisfies the lemma we introduce the notion 
of guiding path. Denoting by T (o) (R^*^) the tangent space of R^*^ at the point z^"'^ , 



^0 

consider the vector 1/ e T (a)(R^'^) pointing towards z with length e{L + A). Condition 
[A3] implies that there exists (a possibly non-uniqiic) injective map a : {1, . . . , 2d} — > 
{ci,...,c„} (where q is the index of the periodic orbits), with a — cr(l?) and non- 
negative coefficients ai, . . . , a2d such that 

2d 

^ = ^£aiJ?,(,)(4'^)). (22) 
1=1 

We refer to the coefficients as guiding times. We denote || the guiding path between 
Zq"^ and z as 

i i 

G{t,zt\z) = + 5^£a,J?.(,)(4")) + £(t -Y,a,)t,^^{zt^) (23) 
i=i j=i 



for £ Xlj=i ttj <t <e Yl'jLi Let us denote by W-^f := |X^(i)(4''0) • • • > ^ ai2d) (4"'') 
the 2(i X 2d matrix whose columns are the 2d Hamiltonian vectors needed to represent 
if as a linear combination of the type (22). Now, condition [A3] implies that 
|^o.(i)(-Zo"'')' • • • ' ^o-(2d)(-2o"'')| is a basis in M^*^, thus is invertible and from (22), 
denoting ~^ = (oi, . . . , a2d), we get 



a\\ < 



V 



< D 



(24) 



At this point we abuse the notation by identifying vector fields and points in 



T>2d 
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where D — sup^^^sup- 



C 



[L + A) < oo depends only on L + A. Hence 



the guiding times are uniformly bounded. 

Without loss of generality ^ we will prove the lemma for d = 1, and we write 
V = eaiXi + £02^2- Then the guiding path consists of two segments. We will show that 
these two segments can be shadowed one at a time, hence the code ^^^^ will be obtained 
by updating the code ^^"^ twice. The rules for updating the code are 



and 



m 



Ck 



Cl 



cib) 
^k 



C2 



for A; < 0, 
for all A; > 0, 



for k < 
for all k > 



ai 



ai 



(25) 



(26) 



where [•] denotes rounding up to the next integer and Tc^ denotes the period of the 
periodic orbit Lc^ (also note that by (24) Oj are of order 0{e)). We define N — 



and A^^ 



and let iV = + NK 



We begin by updating the code according to (25) which gives us the code For 



^b) 



. Let 



any p < we pick the trajectory of (6) corresponding to ^^^^ such that 
us also consider two auxiliary sequences, z^^^ which corresponds to the code ^^^^ and 
z^^ which corresponds to the constant code (ci)°°. Both sequences satisfy the initial 
condition 



iib) 



By the definition (25) we have = Ci for all A; > 0. Then 



Lemma 3 (Lemma 5, [4]) For any Kq > 0,to > 0, there is Eq > such that for any 
\e\ < £o o-nd any two codes and such that for some index j 

to 

^ e _ 

the inequality < eKq implies 



f 1 . = £2 



j+i 



< ^ < No{e) = 



^j+N ^j+N I 



< eCie'^^' 0<N <No{e), 



where 



and 



d(t)c 



d{x,y) 



4r 



1- A 



dz 



d(t)c 



d{x,y) 



max 



dXr, 



dz 



dye 



dz 



% The proof for d > 1 is analogous to the case d=l. The only difference is that the guiding path will 
consist of more segments. Consequently the constant A depends on d. 
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jj+AT - ^|+7v|| < 3eCi 0<N< 



1 



implies that 



< 3s U 



4r 

'■'cici llci I _ \ 



for aRO < k < N. By Lemma 1 we have the following bound 



{&)l 



'0 



'0 ll ^ 



\Zn ' - z^ 'W < eK, 



(27) 



(28) 



since z^^^ and z^^^ share the same code for A; < and z^p'^ = Zp""^ . Now, since i^^^ and 5^^^ 
are generated from the identically same code, Lemma 3 implies that 



<3e (||0cicJlciY^ + ^) 



for alio < k < N. Combining the estimates (27) and (29) we obtain 



'^k ^k 



< 



^b) 



^k ^k 



(b) 



+ 



^k ■^k 



<3£ (^2 ||0,,,J|pi ^ + 



(29) 

(30) 
(31) 



for all < A; < A^. Now it remains to prove that z'^'^ stays close to the guiding path G. 
By definition we have 

G{t, zf,z) = - tQ-VJci(4"^) for < t < eai, 

where Q."^ is the inverse of the symplectic matrix. The map (6) takes a simple form for 
the constant code (ci)°°, see [4] and the proof of Lemma 1 therein. 



(a)^ 



Therefore 



z^~ — G{eai^ z^\ z) <0(£^)+£max sup \Tc{u^v)\\Xc{u,v) 

^ {u,v)ev 

<£(l + max sup \Tc{u,v)\\Xc{u,v)\), 



(32) 
(33) 

where the 0{£^) term is the difference of the two maps and the second term compensates 
for the fact that ai in general is not a multiple of , this round off error is bounded 
by the strength of the vector field times the largest period. Collecting the estimates we 
obtain 



i<| - G{ea,, z^^\z)\ < 3£ [2 ||0c,cJlci + + e 



(34) 
(35) 



+ emax sup \Tc{u^v)\\Xc{u,v)\ , 
{u,v)eD 

which is the accuracy with which we have shadowed the first segment of the guiding 
path. 



Geometric shadowing in slow-fast Hamiltonian systems 



16 



Let us continue by shadowing the second segment of the guiding path. We begin by 
updating the code to using the rule given by equation (26). We pick the trajectory 

which corresponds to ^^^^ and satisfies z^^ = zl"'\ We repeat the arguments above 
to shadow the second segment of the guiding path. Indeed, we consider two auxihary 
sequences zj^^ which corresponds to the code ^^^^ and zl^ which corresponds to the 
constant code (02)°°. Both sequences satisfy the initial condition 

G{eai, z^''\ z) . 

By the definition (25) we have — C2 for all k > N. Then by Lemma 3 

< 3£ ||^C2C2llcl 



Ah) _ cj 
N+k N+k 



A 



(36) 



for all < A; < iV''. Lemma 1 together with (34) yields the following bound 



\z^^ - G{£ai,z^Q\z) 



4r 



< max sup (6£ ||0cc||ci T + 

{u,v)eD 1 - A 

+ £ + £ \Tc{u, v) \ \Xc{u, v)\) . 



Then Lemma 3 implies that 



z^~^ - z^~^ 

N+k N+k 



( 4r 

< 3£max sup 7 ||(/)cc||(^i \ A^K ^\ ^\Tc{u^v)\\Xf.{u^v)\ 

" {u,v)ev\ l-'^ 

for all < A; < N^. Combining the estimates (36) and (37) gives 



(6) _ C2 
N+k N+k 



< 


(b) _ Ab) 


+ 


Ab) _ C2 




N+k N+k 




N+k N+k 



(4r 
8 ||0cc|lci T + + 1 + \Tc{u, v) \ \Xc{u, v) 



(37) 

(38) 

(39) 
(40) 



for all < < A'"^. Using that = G{t^, z^ , z) the guiding path is shadowed by z'~_^^ 



exactly analogously to the first segment. The result is 

2- — G(ea2, zl)^\ z) <£ + £max sup \Tc(u,v)\\Xc(u,v)\ 

{u,v)ev 

Using that -z^-^^^t = z^n ^y definition and combining (39) and (41) gives 

z^N — G{ea2, Zq\ z) = 

( 4r 

<4£max sup I 6 ||0cc||ci ^] r + 3ir + 1 + |Tc(m, v) | |Xc(m, f ) 



(n,i))el> 

Thus, using the definition (11) of A 



•'ccllc'i 



(41) 

(42) 
(43) 



A — 4 max sup ( 6 ||^, 



4r 



A 



+ 3K + 1 + |re(u,t;)| \XXu,v) 
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and the definition of guiding path gives 



< sA. 



The constant A is uniformly bounded on the compact space D. Lastly, from the fact 
that the guiding times ai are uniformly bounded it follows that the length of the guiding 
path is bounded by eC\ where C\ is a uniform constant. This and the estimates above 
implies that 

^^^^ - - < e{Ci + A) iorQ<k<N. 
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